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Edge fracture occurs frequently in non-Newtonian fluids. A similar instability has often been reported at 
the free surface of fluids undergoing shear banding, and leads to expulsion of the sample. In this paper the 
distortion of the free surface of such a shear banding fluid is calculated by balancing the surface tension against 
the second normal stresses induced in the two shear bands, and simultaneously requiring a continuous and 
smooth meniscus. We show that wormlike micelles typically retain meniscus integrity when shear banding, 
but in some cases can lose integrity for a range of average applied shear rates during which one expects shear 
banding. This meniscus fracture would lead to ejection of the sample as the shear banding region is swept 
through. We further show that entangled polymer solutions are expected to display a propensity for fracture, 
because of their much larger second normal stresses. These calculations are consistent with available data in 
the literature. We also estimate the meniscus distortion of a three band configuration, as has been observed 
in some wormlikc micellar solutions in a cone and plate geometry. 



I. INTRODUCTION 

Many complex fluids are dramatically influenced by 
shear flow, which can easily disrupt the slow mi- 
crostructural relaxation of these fluids. Examples in - 
cludc polymer sol utions (Boukan v and Wan 3, 12009ft . 
wormlike micellar dSchmitt et all 11994 ) and lamellar 
( Diat. Roux. and Nalletl. 19931) s u rfacta nt solutions, col- 

and telechelic 
In many cases 



1994) 



2001) 



loidal suspensions ( Chen etal\, 
polymer networks ( Michel et all 
shear flow can induce an apparent transition to a state 
with a different microstructure and apparent viscosity, 
which can lead to macroscopic bands of mat erial that 
coexi s t, much like ph ase separation, in flow ( Olmstedl . 
120081: iFieldinsi l2007t ). In the most commonly observed 
scenario, the system forms two or more 'shear bands'; 
layers of high and low shear rate material (of equal 
shear stress) that coexist at volum e fractions consistent 



with an imposed average shea r rate dSchmitt et al. I.ll994t 
iBritton and Callaghanl Il999t ICappelaere et all \199% . 



As the average shear rate is increased, the width of the 
high shear rate band increases, while a constant stress 
is maintained (in an idealized planar Couette geometry). 
The measured shear stress as a function of applied av- 
erage shear rate is the flow curve, and contains a broad 
stress plateau for shear rates in the banding range. 

Shear banding can result when the constitutive rela- 
tionship between shear stress and shear strain rate, as- 
suming homogeneous flow , has a stress maximum and 

is thus non-monotonic (jSpenlev. Cates. and McLeishl . 

Il993h (see Fig. [St below). Homo geneous flow is un- 
stable for applied shear rates in the region of the con- 
stitutive curve with a negative slope, and this insta- 
bility can be resolved by adopting the shear banding 
state. This constitutive instability is present in the 
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Doi-E dwards (DE) reptation model for entangled poly- 
mers ( Doi and Edwards! . Il989() for sufficie ntly weak levels 
of Co nvected Constraint Release (CCR) ([Graham et all 
I2003D . and w as only recently been observed in poly- 
mer solutions (iTapadia and Wand . 120031: iHu et all 12007 
Boukanv and Wand . 120071 l2009t lAdams and Olmsted 



2009bD . Wormlike micelles have a similar constitutive in- 



stability, due to the comb i nation of reptation and micel- 
lar breakage (ICatesl 19901 Rehage and Hoffmannl . [l99lt 



ISpenlev. Cates. and McLeishl 1993 ). and shear banding 
has b een studied in these systems for decades ( Berretl . 
I2005D . 

In many experiments on shear banding the free 
surface fractures and the sample is ejected from 
the device. Although it is often reported as a nui- 
sance and anecdotally, it is widespread in both 



worm -like mice llar solutions (IBritton and Callaghan , 



1997 



Berret. Porte, and Decruppg. 11997 



Lopez- Gonzalez. Holmes, and Callaghan . 20061) and 



entan gl ed polymer solutio ns ( Inn. Wissbrun. and Dennl . 
120051 : ISui and McKennal . 120071) Fracture and ejec- 
tion can occur at some point on the stress plateau, 
which corresponds to a certain minimum width of 
the high shear rate band. This is eviden t in th e 
sxperiments of iBerret. Porte, and Decruppd ( 1997t ). 



in which they reported surface fracture on the stress 
plateau in a cone-and-plate rheometer. They attributed 
this to a well-known secondary flow instability in 
cone-and-plate rhcometers, due to the balance of the 
second normal stress diffe r ence with su rface tension 
([Tanner and Keentokl . Il98l lLarsonl . fl992l) : 



\N 2 \ > 



2lL 

3W' 



(1) 



where W is the maximum cone-plate separation and 
N% = T yy — T zz is the second normal stress difference. 
The balance of normal stresses with surface tension leads 
to a radius of curvature R ~ 7/-ZV2. If this radius is too 
small, then the interface must curve too tightly to fit 
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inside a wide gap W , and fracture results. 

In this paper we propose a simple model that gener- 
alizes this idea to incorporate shear banding, which can 
also address this lack of surface integrity. We study the 
entire meniscus shape, as determined by second normal 
stresses in the two shear bands and contact angles, and 
find a range of conditions under which the meniscus can- 
not maintain the shape demanded by mechanical equi- 
librium. As the shear bands change size, for increasing 
imposed average shear rate, one or the other band devel- 
ops a width that cannot support the meniscus curvatures 
demanded by the normal stress balances in both shear 
bands and continuity. We establish the conditions for 
mechanically stable bands, and i llustrate this behaviour 
using non- monotonic forms of the ljohnson and Segalmanl 
(I1977D and lGiesekusl (|l982l) constitutive models. 

We compare this explicitly to data in the litera- 
ture. There is significant data on wormlike micelles, 
which are for the most part consistent with our calcu- 
lations. We also comp are with recent work on entan- 
gled polymer solut i ons (Em. Wissbrun. and Dennl . 20051 : 
ISui and McKennal . l2007t iTapadia and Wand . l200l . 
which are now known to shear band and have gen- 
erated much discussion in the literature. We show 
that the larger apparent viscosities of entangled poly- 
mer solutions leads to more unstable shear bands; we 
hope that this helps resolve some of the contradic- 
tory results in the literature as to whether or not 
shear banding is intrinsic to the constitutive behavior 
or due to a meniscus distortion. In our view shear 
banding can initiate the dis tortion and fracture seen 
in some recent experiments dlnn. Wissbrun. and Dennl . 
120051: ISuT and McKennal . I2007D . 

In Section|TT]we present our model, flow geometry, and 
mechanical balance conditions. In Section IIIII we calcu- 
late the meniscus shapes, and derive meniscus integrity 
criteria that depend only on the second normal stresses, 
the gap size, the surface tension, and the contact angle. 
In Section IIVI we compare our calculations with experi- 
ments on wormlikc micelles and polymer solutions, and 
we conclude in Section [V] Appendix [B] collects the rele- 
vant information for the Giesekus and Johnson-Segalman 
models, and Appendix [Cl contains the details for calcula- 
tions in which a center high (or low) shear rate band is 
sandwiched between two low (or high) shear rate bands. 



II. MODEL AND MENISCUS INTEGRITY 
CONDITIONS 

Constitutive Equations — We consider an incompress- 
ible fluid obeying the following relationship between 
shear stress and shear strain rate 



-pi + 2?/D + £ 



(2) 



solvent or other fast modes), and v is the velocity field. 
The stress tensor £ is an additional viscoelastic stress 
that has its own dynamical equation of motion. We will 
illustrate our results using the Johnson-Segalman (JS) 
and Giesekus models, whose details are presented in the 
Appendix. We will consider steady creeping flow, so 



V -T = 



(3) 



where I is the identity tensor, D = | [Vv + (Vv) T ] , p 
is the isotropic pressure determined by incomprcssibility 
(V ■ v = 0), i] is an assumed Newtonian viscosity (due to 



Shear banding usually develops only two bands, 
although occasionally more complex structures are 
seen, such as a three band configuration in cone- 
and-platc (Fig. [T]). The stress gradient in cylindri- 
cal Couette flow typically ensures that two bands de- 
velop with the high shear rate phase near the in- 
ner cylinder ( Olmsted. Radulescu. and Lul . l2000h . while 
the much weaker stress gradient of cone-and-platc 
flow could explain the mo r e com plex structures 
(lAdams. Fielding, and Olmstedl . 120081 ). Alternatively, 
iKumar and Larson! 1)20001 )" showed that unidirectional 
shearing flow, in the cone and plate geometry, of adjacent 
fluids with different normal stresses is incompatible with 
momentum balance. Based on this three band state, we 
calculate the meniscus distortion of two band and three 
band configurations. We specialize to a planar Couette 
geometry for all calculations. We may sometimes refer 
to a particularly shear banding configuration as 'unsta- 
ble' if the mechanical balance condition does not yield 
a physical reasonable interface. However, we empha- 
size throughout that we do not calculate the conditions 
for dynamic instability, but rather the conditions under 
which a meniscus solution is consistent with momentum 
balance. This 'instability' or lack of solution may may, of 
course, be preempted by an instability due to secondary 
flow arising from purely dynamical considerations. 

Shear Flow and Mechanical Balance Conditions- We 
consider steady laminar flow between two parallel plates 
a distance W apart; one plate is stationary and the other 
plate moves with a constant speed V in the x direction 
(see Fig. [5]). The flow gradient and shear rate only vary 
in the y direction, dv x /dy = j{y). The applied average 
shear rate is 7 app = V/W . The free surface (meniscus), 
specified by a height h(y), is in the z direction, and we 
assume no-slip boundary conditions at the solid walls. 

If the free surface is considered at all it is usual for the 
meniscus to be flat (parallel to the plane z = 0). Cur- 
vature of the free surface is problematic; stress balance 
at the surface implies that the base flow near the sur- 
face will, generally, no longer be given by v = (72/, 0, 0). 
For h! <C 1 this may not be a problem; otherwise sec- 
ondary flows will develop and, generally, we expect hy- 
drodynamic instabilities to preempt our estimates of non- 
existence of surface integrity. This is discussed in the 
Appendix. 

The fluid is assumed to have a non-monotonic constitu- 
tive relation, with shear bands that form at shear rates 71 
and 72 two shear bands Fig.[3ta). The shear band widths 
wi and w 2 are given by w 1 = -(7 app - 72)/(72 _ ^ and 
u>2 = (7a PP - 7i)/ (72 " 7i) j where here and through- 
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FIG. 1. (Left) Flo w birefringence image of shear banding in a cylindrical Couette rheometer; the bright band is flowing at the 
higher shear rate (|Cappelaere et al. Ul997l ). (Righ t) NMR image of shear bandi ng in a cone-and-plate rheometer, with a high 
shear rate central band (white or black) ( Fig. 2a lBritton and Callaehanl (jl997l )). A free exterior fluid surface is specified. 



out this paper all lengths with a carat " been normalized 
relative to the plate gap size W. The shear bands are 
assumed to partition, and vary in thickness, along the 
flow gradient direction. 

Stress balance (Eq. implies that T xy and T yy are the 

same in each shear band. This implies that = E^, 
where ^ and ( 2 ) refer to the two shear bands. Moreover, 
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Since the normal stresses will generally be different in the 
two bands, the pressures pi and p 2 will differ. These pres- 
sures must then balance, together with E zz , against the 
curvature of the meniscus and the atmospheric pressure 

Patm • 



length for the shear banding configuration: 
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The three equations Q E]) relate four unknown quan- 
tities: the pressures and meniscus curvatures in each 
band. By eliminating the pressures we can relate curva- 
ture radii, but more information is needed to absolutely 
determine the shape. This will follow by constructing 
a continuous and smooth meniscus. The balance at the 
meniscus, Eq. (JQ, determines the pressure in each band 
in terms of the meniscus curvature. We will find below 
that the curvatures must change in order to maintain a 
physical meniscus, which thus determines the pressure in 
each band. 



III. MENISCUS SHAPE AND INTEGRITY 



Patm 
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where Ri is the radius of curvature of the meniscus in the 
zth band and 7 S is the surface tension. From Eq. (UJ) the 
difference in T zz between the two bands is given by the 
difference in the second normal stress differences, TiV — 



T S = ^2 - ^2 • Making use of this and Eq. ©, 
we can relate the second normal stress differences to the 
radii of curvature of the two bands 
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where A7V 2 = 7V 2 (2) - N^> . This can be easily calculated 
for a given constitutive relation, and together with the 
surface tension defines a characteristic 'elasto-capillary' 



r (i) 



A. General Shape (two bands) 

We ignore deviations due to complex flows near the 
contact line, and approximate the meniscus of each band 
as the arc of a circle of radius i?,, which may be posi- 
tive or negative. We demand continuity of the surface 
and its tangent at the interface between two bands. The 
conditions above can be fulfilled either by: 

(i) the contact angles being fixed and the fluid adopt- 
ing a height H difference between the contact lines 
along each plate ( movement of the contact l ines h as 
been observed by I Crawley and Graesslevl ( 19771 )): 
this height difference will vary in response to the 
widths of the shear bands; 

(ii) the contact lines being pinned such that H is fixed 
and the contact angles vary in response. 

In both cases the meniscus profiles are governed by the 
same equations. 
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FIG. 2. (Left) Schematic diagram of planar Couette flow. The plate in the y = plane is stationary whilst a parallel plate a 
distance W away moves with velocity V . The free surface at large z is open to the atmosphere. (Right) Velocity profiles for an 
applied average shear rate 7 app = (a) unhanded or (b) shear banded. For banded flow 71 = ^- and 72 = v ~ v , where wi 
and W2 are the widths of the two shear bands. 




FIG. 3. (a) Flow curve (thick line) and constitutive curve (thin line) for shear banding flows. Banding occurs on the stress 
plateau at T* y , for applied shear rates 7 app such that 71 < 7 app < 72. A 'lever' rule relates the widths of the shear bands to 
the applied shear rate, wi = 1^(72 — 7a PP ) / (72 —71) and W2 — W — wi . (b) Profile of the fluid surface (meniscus) between flat 
plates (flow in the x direction). Each band has a circular profile, and the surface is continuous and differentiable; <j>\ and 02 
are the contact angles at the two plates. The contact lines between fluid and plates move up and down as the widths of the 
shear bands alter; H is the height difference between the contact lines. Note that the curvatures may be positive or negative. 



The height profiles of the two bands are given by 



Kv) = 



Ri \ \ 1 1 - l jj- + 




y-W 



(0 < y < vii) 



(8) 



cos <j>2 I - sin (j> 2 + H {wi < y < W) , 
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where 0i and 02 are the contact angles at either wall, and h(0) = and h(W) = H . Continuity of h at y = wi 
requires 



H = (Ri - R2)\j 1 — ( -j^- + cos 0i ) + i?2 sin 2 - R\ sin 0i , 



and continuity of h! at y — w\ requires 



wi | w 2 
Ri R-2 



COS 01 — COS 02 



(9) 



(10) 



B. Fixed Contact Angles 



We allow the surface tensions at the fluid, walls, and atmosphere interface to determine the contact angles and 
assume that the contact angles persist through the onset of banding. Then Eq. ((SJ) and Eq. (jTU)) completely specify 
the shape: 



Ri 



1 



- COS 01 — COS 02 — W-2.A 



and R 2 = 



1 



- COS 01 — COS 02 + W\A 



where the dimensionless distortion parameter 



A = 
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(ii) 



(12) 



controls the shape. In the limit of high surface tension, \A\ ~ 0, both radii are equal and completely determined by 

the contact angles. We have chosen negative values for A here since we expect N 2 < < in the high shear rate 
band for most polymer and micellar solutions (a similar analysis can be done for A > 0). In Section HVl we analyze 
recent experiments and estimate —A ~ 0.8 - 3 for shear banding wormlike micelles and —A ~ 3 - 140 for entangled 
polymer solutions. 

I 



C. Integrity of Meniscus 

An interface solution exists under two conditions, one 
mathematical and the other practical: 

• The curves should not have infinite slope or pass 
through each other. This leads to the condition 



W2 COS 01 — Wl COS 02 — Wl W2 A 



< 1 



(13) 



• The increase H in height across the gap should not 
be large enough for the sample to climb out of the 
cell. This obviously depends on the loading condi- 
tions, and can generally preempt the mathematical 
condition above. 



D. Equal Contact Angles 

We first consider equal contact angles 0i 
for which the meniscus integrity condition is 



\(W2 — Wl) COS0 — WiW2A\ < 1, 

where the criterion 



1 — cos < — < 1 — COS ( 
Ri 



(14) 
(15) 



is required for Eq. © to yield real solutions. Fig. @Ji 
shows regions of integrity for A = —3.5 and two cho- 
sen cases. For a wall and fluid combination that fixes 
= 60° the entire shear stress plateau is accessible and 
supports a integral meniscus. However, for a wall and 
fluid combination that fixes = 40° some values of W2 
are not allowed and the corresponding applied shear rates 
are not accessible on the stress plateau. For small con- 
tact angles the low shear rate portion of the plateau is 
accessible, while for high contact angles the high shear 
rate portion of the plateau is accessible. Fig. 0b shows a 
possible flow curve for the = 40° case. 

To understand the lack of integrity of the meniscus 
we must examine the meniscus shapes. These are shown 
in Fig. [5] as a function of increasing the applied shear 
rate 7 app across the stress plateau, for both a integral 
case (0 = 60°) and an non-integral case (0 = 40°), for 
A = —3.5. Note that banding first initiates when the high 
shear rate band develops a non-zero width, i.e. when W2 
grows from zero for 7 app > 71 . 

Consider first the = 60° case (Fig. |5k). At w 2 = 0.2 
the contact line of the higher rate shear band has dropped 
while that of the lower shear rate band has risen. Both 
surfaces have negative radii of curvature. By W2 = 0.4 
the radius of curvature of the lower shear rate band has 
become positive. The height difference between the con- 
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FIG. 4. (a) Integrity contours for equal contact angles cf>2 = <j>x = (f>, for A = —3.5. For a wall and fluid combination that fixes 
4> = 60° all shear band widths are allowed and the applied shear rate will traverse the plateau. For a wall and fluid combination 
that fixes (f> — 40° not all shear rates accessible; the surface of the fluid becomes fracture for some widths W2 of the high shear 
rate band satisfying < W2 < w^- All widths allow meniscus integrity for ] cos</>] < cos0* = 0.661. (b) Flow curve (thick 
line) for a wall and fluid combination with cf> = 40° , showing the inaccessible segment of the stress plateau. 



tact lines increases and reaches a maximum at w 2 — 0.65, 
but the surface still maintains its integrity. For larger 
w 2 the height difference decreases until shear banding 
ceases at w 2 — 1 (7 apP = 72)- For the <f> = 40° case, 
however (Fig. [5b), the integrity of the surface cannot be 
fulfilled when the width of the high shear rate band sat- 
isfies 0.6 < w 2 < 0.83 (see Fig. [SJi).As can be see in 
Fig. [5t>, the interface develops a kink immediately be- 
fore the onset of fracture. In the region of fracture we 
find that there is no continuous solution to the meniscus 
profiles. 

Figures [5]^cd) show the height difference H, and the 
surface curvatures for the different shear bands ^ nor- 
malized by the respective band size, as a function of shear 
band width w 2 (or equivalently the applied shear rate 
Tapp)- The surface will maintain integrity so long as so- 
lutions for h'(y) (and consequently H) exist. The dotted 
lines indicate between which values 7^ and must lie 

Hi R2 

in order to maintain integrity. 

The integrity contours for a range of distortion param- 
eters A are shown in Fig. [51 and the different regions 
of integrity are summarized in Table fl] For A > — 2 
the meniscus is always integral. For —4 < A < —2 the 
meniscus is integral for contact angles satisfying 



cos 2 6 < cos 2 1 



= \A\ 1- 



(16) 



For large magnitude \A\ > 4 the meniscus looses in- 
tegrity at all contact angles, for some regions of the stress 
plateau. For A < 0, as in our case, this occurs when the 



width of the high low shear rate band is between the 
limits 



1 



2cos< 



i± 4/1- 



A\A\ (1 + cos< 
(\A\ + 2cos^ 



(17) 

These limits also apply for —4 < A < —2, for | cos(f>\ > 
I cos 6* I . 



E. Different Contact angles 

For different contact angles the situation is more com- 
plex. For A < the meniscus remains integral for all w 2 
for A > -8 and 



COS 02 > 



V| COS 0! - 1| - y/\A\ -1. 



For A < — 8 the meniscus fractures somewhere along 
the stress plateau for any combination of contact angles. 
This condition is illustrated in Fig. which shows the 
regions of integrity and fracture as a function of both 
contact angles, for A = —3.5 (Fig. 0k is a slice through 
this figure for <p\ = cf> 2 ). Fig. [TJd shows that an asym- 
metry in contact angle increases the region for meniscus 
fracture. For <j>i = 0, corresponding to a wetting surface, 
a larger contact angle leads to a fractured meniscus for 
I .A I < 2, which would always be integral for equal contact 
angles. 
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Meniscus profiles for ^ = 60° 



Meniscus profiles for <j) = 40° 
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Height difference H (w 2 ) for <f> = 60° Height difference H(w 2 ) for (j) = 40° 





FIG. 5. (a,c) Meniscus profiles h(y); and (b,d) height difference H and curvatures normalized by the gap size, ^, as a function 
of the high shear band width 1S2; for A = —3.5. The dotted lines are the integrity limits on ^ from Eq. (|15[) . For the <f> = 60° 
case (a,c) the surface remains integral for all shear rates in the plateau region, while for the (f> = 40° case (bd) the surface of 
the fluid will fracture on some portion of the stress plateau. 



F. Pinned Contact Lines R%, R2, 4>i and <p2] the fourth condition is the require- 

ment of volume conservation for an incompressible fluid. 
For pinned contact lines we ascribe a particular value Fi §- IE is a typical example in which the contact lines 
to H . Eqs ©, © and (|10[) arc not sufficient to specify 
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FIG. 6. Regions of integrity for equal contact angles cj> and for different values of the distortion parameter A. 



Range of A Integrity of meniscus 

A < — 4 Fractures for all contact angles for < u>2 < ■ 

—4 < A < —2 Fractures for |cos<^>| > |cos0*| and < W2 < ■ 

—2 < A < Integral for all contact angles <j>. 

TABLE I. Criteria for (lack of) integrity as a function of the distortion parameter A, for equal contact angles c/>i = cj>2 = </>■ 



have been pinned so that their difference in height H is 
zero. For applied shear rates associated with the stress 
plateau the meniscus contorts but retains integrity for 
all shear rates when —8 < A < 0; for A < — 8 not all 
shear rates are accessible and the meniscus will fracture 
at some shear rate. 



G. Three Bands 

Fig. [T] shows an example of three shear bands, 
visualized in a c o ne-an d-plate rheometer by 
iBritton and Callaghanl ( 1997t ) using NMR velocime- 
try. Three bands have not been observed in cylindrical 
Co uette flow, and this differen c e wa s rationalized 
by lAdams. Fielding, and Olmstedl ( 2008| ) as due to a 
combination of boundary conditions that favor the low 
shear rate phase and the relatively weak stress gradient 
of cone-and-plate flow. Motivated by this result, we 
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FIG. 7. (ab) Regions of integrity as a function of contact angles 0i and cj>2, and the width iZ>2 of the high shear rate band for 
distortion parameter A = —3.5. The meniscus is integral for any combination of W2, </>i, and <j>2 in the integral region. Fig. [4^. 
is the intersection of this integrity surface with the plane 4>i — <j}2- For A = — 1 and <f>2 < 90° the meniscus is integral for all u>2 
(i.e. across the entire stress plateau). For A < A* = — 2 the meniscus fractures at some shear band width W2 for all contact 
angles <f>2. (c) Integrity contours for 02 = 90°. The meniscus fractures at some W2 for all 02 if A < —(3 + 2\/2). (d) Contours 
specified by A* as a function of contact angles (4>i,4>2) that enclose regions of integrity across the whole shear stress plateau. 
On the dashed line cj>i = 02- For A* = —3.5 the dotted arrows indicates the contact angles for which the meniscus will be 
integral for all W2- The dotted lines indicate fracture at some W2 as shown in Fig. [4^,. 



have calculated the distortion of the meniscus in such 
a configuration and find that a three band state does 
indeed allow meniscus integrity, for certain ranges of 
parameters. This analysis is given in Appendix [Cj 



IV. COMPARISON WITH THE LITERATURE 

A. Theoretical Method 

We examine existing data in the literature, which ex- 
hibits both stable shear banding across the entire stress 
plateau and an instability such that the entire stress 
plateau could not be accessed; we then assess whether 
or not meniscus fracture is expected, based on an esti- 
mation of the distortion parameter A — W AN2/j s - We 
fit the data from the flow curves, including the stress 
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FIG. 8. Stability regions of profiles for contact lines pinned such that H — 0. (a) The full width of the shear stress plateau is 
accessible for —8 < A < 0. For A < —8 the surface will fracture at some applied shear rate < 7*. (b) Meniscus profiles for 
A — —8. The surface maintains integrity for applied shear rates up to 7* at which rate the high shear rate band has achieved 
u>2 = 0.3 of the width of the gap. 



plateau and the shear rates 71 and 72 in the two shear 
bands, to both the Giesekus and Johnson-Segalman (JS) 
constitutive models. From the models we can evaluate 
the predicted second normal stress difference N2 in the 
two shear bands (see Appendix |B} , while the gap size 
W is taken from the experimental conditions and the 
surface tension 7^ is estimated using literature values. 
With this in hand we can then compare the stability or 
instability of banding states, inferred from the experi- 
ments, with our calculations. Since the contact angles 
are unknown, we can only determine whether our crite- 
ria for instability are consistent with the numerical values 
for A (assuming equal contact angles). We are also lim- 
ited by the quality of the available constitutive models: 
neither the JS or Giesekus models are expected specifi- 
cally apply to wormlike micelles, but they can support 
shear banding and have non-zero second normal stresses. 
These models, developed for polymer melts, may be 
more applicable to semi-dilute poly meric solutions, which 
have a n explicit Newtonian solvent. Fischer and Rehagej 



(Il997l). lYesilata. Clasen. and McKinlevI ([20061 ) . and 
iHelgeson et al. ( 20091) show that the steady state non- 
linear rheology and shear thinning of worm-like micelles 
is well described by the Giesekus model; we offer the 
Johnson-Segalman model for comparison. In fact, the 
two models yield very similar quantitative predictions 
and thus do not substantially differ insofar as determin- 



ing the properties of the meniscus. 



B. Wormlike Micellar Solutions 

Wormlike micellar solutions of both cetyltrimethylam- 
monium (CTAB) and cetylpyridinium choloride/sodium 
salicilate (CPCl/NaSal) have been extensively studied. 
Fig. IH1 s hows the data o f a C TAB solution, as mea- 
sured by IHelgeson et all (|2009l) and fit by them to the 
diffuse Giesekus model. In this case the entire stress 
plateau was accessible. We have also fitted it to the non- 
diffuse Giesekus and JS models. Strictly, one should fit 
the stress plateau using a non-local (or diffuse) model 
(|Lu. Olmsted, and Balll l2000t ); however, our fits ob- 



tained by choosing the stress plateau 'by hand' differ 
insignificantly from more precise fitting. Hence we use 
a local model for the remaining fits in this paper. The 
Giesekus model fits the high shear rate branch much bet- 
ter than does the JS model. 

The fitting parameters are shown in Table HU The 
difference between the diffuse and non-diffuse Giesekus 
models values for AN2 is less than 3%. We estimate the 
distortion parameter to be A ~ —0.8, which is well within 
range (A > —2, from Table A} for which we expect to find 
stable shear bands across the entire stress plateau for all 
contact angles. 
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FIG. 9. (a) Flow curve fo r a CTAB solut i on fit ted to the (numerical) diffuse Giesekus model, measured in a Couette rheometer 
with a gap W = 0.5 mm (|Helgeson et 01.1,120091 ). Fits to the non-diffuse Giesekus (b) and JS (c) models (fitting parameters are 
given in Table llTjl . 
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T/°C 


f/mm 


G/Pa 


r/s 


a e 


ii/s- 1 


72/s -1 
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n 


A 
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0.5 
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65 
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0.43 0.022 


40 
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0.5 
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0.11 


0.88 0.015 


11 


60 
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-0.9 Gi 
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0.09 


0.41 0.056 


8 


90 


-0.48 


-1.2 JS 


Fig. OH 


CPC1 


240 


20 


0.5 


260 


0.94 


0.87 0.002 


1 


54* 


-0.60 


-2.2 Gi 




12% 








268 


0.88 


0.38 0.010 


0.85 


54* 


-0.49 


-1.9 JS 



TABLE II. Parameters used to fit the flow curves in Figs. l9l ITOl and llll for the Gies e kus ( Gi) or Johnson-Segalman (JS) models. 
We estimate a surface tensi on 7 S ~ 37mN/m for aqueous CTAB (| Christian et all 1 19981 ) and 7, ~ 32 mN/m for CPCl/NaSal 
<|Akers and Belmontd . [200^ V For the stable solutions (CTAB and 8% CPC1) the shear rate 72 of the high shear rate band was 
ta ken from the end of the stress platea u, while for the unstable solution (12% CPC1) 72 was estimated from the scaling in Fig. 8 
of lBerret. Porte, and Decruppd (| 19971 ). 



We have found a few experiments that show 
a clear instability as the stress plateau is tra- 
versed, or that cannot reach the end of the stress 
latea u. B erret and co-workers rtBerret. Roux. and Portd 
19941 ) and lBerret. Porte, and Decruppd (|l997l )) studied 
CPCl/NaSal solutions at different concentrations and 
temperatures, in a cone-and-plate rheometer with a free 
surface. They found stable stress plateaus at micellar 
solutions close to a non-equilibrium critical point, where 
the difference between the shear banding phases vanishes. 
Farther from the critical point, where the shear banding 
phases are more distinct, and hence one expects AN2 to 
b e larger, they report unstable stress plateaus (see Fig. 6 
of lBerret. Porte, and Decruppd (|l997t )). Fig. ITU1 show the 
data and fits for an 8% solution at T = 30° C, which was 
close to the non-equilibrium critical point. The fitting 
leads to a distortion parameter between A ~ —1.0 and 
A ~ —1.3, which is consistent with the meniscus main- 
taining integrity across the entire stress plateau. 

The data for an unstable CPC1 solution are shown 
in Fig. [TTJ at 12% and T = 20° C. In this case 
the stress plateau can only be traversed as far as 
7a PP = 7.1s" 1 , at which point the fluid became un- 
stable ( Berret. Porte, and Decruppel Il997f ). To calcu- 



late the distortion parameter we require the shear rate 
72 in the high shear rate phase, which we estim ate 
from Fig. 8 of ( Berret. Porte, and Decruppel . Il997l ) to 
be 72 — 54 s -1 . This is consistent with the mea- 
surements of the high shear rate branch perfo rmed by 
lLopez-Gonzalez. Holmes, and Callaghanl ( 20061 ) on the 
same material with a free surface: Fig. 11 of their paper 
shows a a 10% sample at 25° C. Fits to the Giesekus and 
JS models respecti vely yield A ~ —2.4 and A ~ —2.0, fo r 
the 12% sample of lBerret. Porte, and Decruppd ( 19971 ). 
at T = 20° C. In the range A < —2 we expect menis- 
cus fracture for some contact angles, and across some 
region of the stress plateau. The experiments show an 
instability at W2 = 0.12. Fig. [5] shows that such a small 
value of W2 = 0.12 and A < — 2 would be consistent 
with a contact angle near 180°, or complete wetting. 
We caution, however, against using the numerical results 
from applying these fairly crude constitutive models. It 
is clear, however, that the distortion parameter predicts 
that the 12% fluid at T = 20° C should be more unstable 
than the 8% fluid at T = 30° C. We have also esti- 
mated the dist ortion parameter from the closely related 
experi ments of lLopez-Gonzalez. Holmes, and Callaghanl 
( 20061 ) (see Fig. 5 from that paper), for which we esti- 
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mate —A ~ 2.3 — 2.6, again just into the predicted range 
of possible meniscus fracture. 

We have collected much of the available data in the 
literature in Table Mil In many cases the full banding 
profile was found (noted as 'stable' in Table HIT)) . and in 
all of these cases we calculated A to be in the range A > 
—2, the meniscus retains integrity (Table H}. Similarly, 
all of the unstable data that we have found corresponds 
to A < —2, for which meniscus fracture is expected for 
some contact angles (Table fl}. 



C. Polymer Solutions 

Wang and co-workers revived the experimental study 
of entangled polymers wit h new experiment s that 
clearly sh o w she ar banding (ITapadia and Wand . 12005 



Hu et all l2007t iBoukanv and Wand . 120071 1200 



Tapadia and Wand 's first experiments, on polybutadiene 



(PBD) entangled in its own oligomer, suggested shear 
banding by what they referred to as an 'entanglement dis- 
entanglement transition' (EDT); this data is also broadly 
consist ent with theories based on th e Doi-Edwards tube 
model (|Adams and Olmsledl . l2009bllal ; IWand. l2009t) . In 
the earliest experiments they apparently used a free sur- 
face, and commented that some small edge fracture or 
instability could have occurred, but that its effects were 
negligible. 

While attempting to r eprod uce these experiments, 
linn. Wissbrun. and Dennl (120051) found significant edge 
fracture and instability in the region of the transition 
and the stress plateau, and concluded that the surface 
instability and associated mass loss was responsible for 
the ef fects. This was then studied by ISui and McKennal 
(120071) using cone-and-plate rheometers with gap sizes 
at the rim of W = 0.267 mm and W = 1.363 mm. 
They found significant edge fracture or instability ac- 
companying the shear banding transition, but their re- 
sults suggested that the surface instability could be a 
consequence, rather than a cause, of the shear banding 
tr ansition (or EDT). This is seen most clearly in Fig. 5 
of[I ui and McKennal (j2007l) . which shows that the stress 
plateau can be traversed when a plastic film is used to 
suppress the instability at the surface (Philips and Wang 
perfor med the film-suppressed experiments). ISchweizer] 
(12007D also found that surface deformation accompanied 
the banding, or EDT, in s i milar mat erials. More recently , 
iRavindranath and Wane] (|2008l ) and lLi and Wand (|2010l ) 
have verified that the shear banding (or EDT) transi- 
tion can persist when even when surface effects arc sup- 
pressed. 

To address these issues we have est imated the distor- 
tion p arameter for the experiments of ISui and McKennal 
( 20071 ) on polybutadiene solutions, and on solutions of 
polyisobutylene in pristane (SRM 2490) (Table HVJ). The 
distortion parameter for PBD is quite large and nega- 
tive, between A — —120 and A = —20 for these experi- 
ments, such that instability should be present early into 



the stress plateau for all contact angles. This is consis- 
tent with the experimental results noted above: shear 
banding can be observed if the meniscus is sealed, and 
an instability occurs otherwise. 

On the other hand, we estimate the lower viscosity ma- 
terial SRM2490 to have a muc h smaller A. These value s 
are consistent with the data o f S ui and McKennal (I2007D . 
who noted explicitly that reducing the rim gap delays the 
onset of instability and expulsion as does reducing the 
modulus (G); both of these changes reduce \A\. They 
report controlled strain rate data for the larger 6° cone 
angle (and wider gap), for which we estimate A ~ —30, 
which is well within the regime where one expects frac- 
ture; see Fig. [12j They observe noticeable meniscus dis- 
tortion at 7 ~ 20 s -1 , and fluid expulsion at 7 ~ 50 s -1 . 
This suggests only a very small window of stability on 
the stress plateau, which is consistent with our calcula- 
tions. For the smaller 1° cone angle they only perform 
controlled stress experiments, and did not observe mass 
loss at the accessible strain rates (~ 95 s _1 ). We estimate 
values of A ~ -4.9 (JS) or A = -6.2 (Gicsekus), which 
is just beyond the limits where we predict meniscus in- 
tegrity across the shear plateau. Hence, their controlled 
stress experiments can access a large part of the stress 
plateau; moreover, we expect that an fracture would oc- 
cur before the end of the platea u is reached, based on 
earlier work on the same fluid by ISchultheisz and Leighl 
(120021) . In their discussion they noted: 



"The conditions at the edge of the cone and 
plate [gap 0.49 mm] can impact the measure- 
ments in several ways, but these effects are 
not easily quantifiable. Perhaps the most sig- 
nificant difficulty is that the Quid can escape 
from between the cone and plate. One indica- 
tor of loss of fluid would be a decrease in the 
moment with increasing shear rate. This de- 
crease was only observed in three experiments 
at 0°C in the step from a shear rate of 63 s _1 
to a shear rate of 100 s _1 . Those three mea- 
surements were discarded. The only other 
evidence of edge effects occurs at the three 
highest shear rates at all three temperatures 
[0, 25 and 50° C], where there is an increase in 
the relative scatter of the viscosity data. For 
this reason, data at the three highest shear 
rates [40, 63 and 100 s _1 ] arc provided as ref- 
erence data only, since the sample geometry 
might not match our assumptions, and the 
uncertainty in t he data cannot be completely 
quantified." — i Schultheisz and Leigh , \200$l, 
p. 22) 

We estimate the stre ss plateau to end at 7 ~ 30 s _1 . 
The measurements of ISchultheisz and Leighl (|2002h . with 
a cone and plate gap of 0.49 mm, thus show that meniscus 
distortion has starte d by 40 s -1 . Even a llowin g for the 
narrower gap of the ISui and McKennal (|2007l ) 1° cone 
and plate, we would expect their fluid to fracture before 
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FIG. 10. Data and fitted constitutive curves for the experiments of iBerret. Porte, and Decruppd (|l997h on CPC1 at 8% (by 
weight) and T = 30° C, using the (a) Giesekus (b) JS models. The geometry was a cone and plate rheometer. Parameters are 
shown in Table PHI 
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FIG. 11. Data and fitted constitutive curves for the experiments of lBerret. Porte, and Decruppd [jl997h on CPC1 at 12% (by 
weight) and T = 20° C, using the (a) Giesekus (b) JS models. The geometry was a cone and plate rheometer. Parameters 
are shown in Table HI The fluid became unstable at the rightmost data point ■, at -yapp = 7.2 s . The open circle o 
shows our estimate of the shear ra te in the high shear rate branch based on extrapolation from the rheology of Fig. 8 of 
IBerret. Porte, and Decruppd (1 1997 ) , tog ether with the temperature dependences in Figs. 5 and 7 of the same reference and 
Fig. 4 of lBerret. Roux. and Portd l| 19941 ) . 



achieving a shear rate of 300 s 1 . 



V. CONCLUSION 

A. Discussion 

Our meniscus distorti on calculation leads to a result 
that is similar to that o f lKeentok and Xud (|l999t) . which 
was not devised for shear banding fluids. They consid- 
ered a more detailed calculation that incorporated the 
flow field due to a small perturbation, but also led to a 
stability condition that balanced the normal stresses with 



surface tension. By contrast, we have effectively assumed 
that the radius of curvature is large enough, compared 
to the gap size, so that one can focus entirely on the 
integrity of the free surface. Hence, we expect our con- 
dition for fracture of the meniscus to be preempted by 
secondary flows as the radius of curvature decreases. 

It would be difficult to experimentally probe the depen- 
dence on contact angle, because entirely different sam- 
ple cells or surface preparations would be needed. How- 
ever, for highly viscous materials the contact angles can 
be se t by the loa d ing p rotocol ( Schweizer and Stocklil . 
120081 ; ISchweizeil 120071 ). and experiments performed 
much more quickly than the true equilibrium contact 
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Oappolanro, Rorrot, Den-urine Oossclv, and Lindner (1997) 



Bcrrc t, Porte, and Decruppe! 1 19971 V The shear rate in the high shear rate branch for the 12% solution was estimated by extrapolating 



from Fig. 8 in this reference. 
Hclgcson. Vasauez. Kalcr. and Wagner (2009). 



Salmon, Colin, Manneville. and Moli 



Lopez-Gonzalez. Holmes, and Callaghan 



'2003) 
Ill2006f). 



eytrar.nlatir.or frnm Fior S in herret Porte and TWrnru^ I 1997f) 



The shear rate in the high shear rate branch for the 10% solution was estimated by 



Lerouge. Fardin. Argentina. Gregoire. and Cardoso ! 20081 ) . Although the data in this reference are apparently stable because the 
entire stress plateau is traversed, the cell top was covered: for an open cell with a free surface an instability occurs at 7a PP — 40 s 
(Lerouge, private communication). 

TABLE III. Values of the distortion parameter A calculated from experimental studies on worm like micellar solutions using 
either the Giesekus (Gi) or JS models. "Stable" experiments accessed the full stress plateau, while in "unstable" cases an 
instability occurred at the apparent shear rate j^pp for a high shear rate band of width w?. The first column identifies th e 
figure above that demonstrates the fit, and the weight fraction use d for the experiments o f lBerret. Porte, and Decruppd l| 19971 ). 
We estimate a surface tensi on 7 S ~ 37mN/m for aqueous CTAB (|Christian et all Il998l ) and 7s ~ 32mN/m for CPCl/NaSal 
(jAkers and Belmontd . I2OO6I ) . 



angle can be reached. Hence, one could systematically 
vary the contact angle to qualitatively test our predic- 
tions (e.g. Fig. [7]). Other more obvious experimental 
tests would be to change the gap size and directly ob- 
serve the deformation of the meniscus as shear banding 
proceeds. 



Although our calculation was performed with shear 
banding fluids in mind, our predictions could be tested on 
simpler fluids. For example, two immiscible fluids with 
different second normal stress behaviors could be pre- 
pared as discs of different thicknesses below T g , loaded 
into a cone and plate rheometer, and then brought into 
the melt state before shearing and observing the menis- 
cus. 



B. Summary 

Motivated by edge-fracture-like instabilities that occur 
during shear banding, we have calculated the distortion 
of the free meniscus in the gradient shear banding con- 
figuration. The second normal stress difference between 
two shear bands determines the radii of curvature of each 
band, and the conditions of continuity and smoothness 
of the interface lead to integrity conditions for the in- 
terface. The integrity limit is defined to occur when the 
meniscus is entirely vertical and thus on the point of over- 
hanging itself; this is probably a conservative estimate of 
meniscus integrity, since we have not explicitly calculate 
the detailed velocity field in the region of the surface. 
It must be stressed that we have presented a simplified 
calculation ignoring the complicated flow field and shear 
conditions near the surface. The necessary calculation 
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FIG. 12. Data and flow curves for the experiments of lSui and McKennal (|2007h on SRM 2400, using a 6° cone angle. ■ indicates 
rate-controlled and 9 stress-controlled rheometry. Fits are to the (a) Giesekus and (b) Johnson Segalman models. 




FIG. 13. Data and flow curves for the experiments of lSui and McKennal (|2007r i on PBD under controlled strain rate conditions, 
fitted to either the Giesekus (a) or JS (b) models. The Giesekus fits show several different fits, to illustrate the difficulty 
of precise characterization of the constitutive models. Table IIVI shows that the distortion parameter A is only very weakly 
sensitive to these different fits. Here, A indicates 1° gap and • indicates 6°. 



has been outlined in Appendix |A"1 

We have calculated integrity diagrams as a function 
of contact angle and the average applied shear rate, for 
given values of the distortion parameter A = WA^/js- 
Since the sample will climb up one wall to to maintain 
integrity, expulsion would usually be expected to occur 
before our theoretical integrity limit has been reached. 
The integrity diagrams depend only on the values of A 
and the contact angles, regardless of the constitutive be- 
haviour. 

By comparison with parameters extracted from exper- 
iments, we find that wormlikc micellar solutions are often 
expected to have stable shear bands, although they can 
sometimes be unst able; this is consistent with t he ex - 
perimental record. iBerret. Porte, and Decruppd (Il997l ) 



found that micellar solutions became more unstable when 
the shear bands were more different from each other 
(i.e. farther from the non-equilibrium critical point in 
concentration-temperature space), and at lower temper- 
atures where the modulus, and hence and thus A, are 
higher. 

Entangled polymer solutions, on the other hand, are 
predicted to easily become unstable due to their much 
higher normal stresses. This is consistent with the re- 
cent body of experiments that show shear banding in 
entangled polymers: unless the surfaces are protected or 
the meniscus shielded, the meniscus d evelops an instabil- 
ity (linn. Wissbrun. and Demi 120051: 1 ui and McKennal . 
I2007T) . 

We also considered the meniscus integrity of three 
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0.68 


0.001 
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50 
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JS 





TABLE IV. Parameters used to fit the flow curves in Figure 112b , and | 12b to the Johnson -Segalman (JS) and Giesekus 
(Gi) and evaluate the distortion parameter A for the study of ISui and McKennal l |2007t ) of polybutadiene-in-oligomer 
(PBD) and polyisobutylene-in-pristane (S RM 2490) solutions. The surface tensi on of pristane is around 26 mN/m 
(|Bascom. Cottington. and Singleterrvl 1964) and that of 1,4 PBD is around 30mN/m (|Brandrup and ImmerguU I1989T ). The 
PBD was fit to the Giesekus model for several possibilities to illustrate the robustness of the parameter A. *For the SRM 
solutions the shear rate 72 of the high shear rate band was estimated based on the parametrization of the shear thinning 
behaviour of the low shear rate branch and the stress plateau. The parameter A is very insensitive to the precise value of 72: 
values of 72 from 200 to 600 s _1 change A by less than 10%. 



bands, as was seen by iBritton and Callaghanl (jl999l) 
in a cone and plate rheometer. [A cone and 
plate is a good candidate for three bands because 
of it s the very weak stress gradient du e to curva- 
ture (jAdams. Fielding, and Olmstedl . I2008T) .] To resolve 
whether or not the three bands can be symmetrically ar- 
ranged, we calculated the surface energy (proportional to 
the surface tension) . The stability of this energy leads to 
a rich variety of possible configurations: the high shear 
rate band can be stable in the center, off-center, or the 
three band state is unstable. 

Since the second normal stress difference is difficult to 
measure, we used the Johnson-Segalman and Giesekus 
models to evaluate it for the experimental examples of 
micellar solutions (low shear modulus) and polymer so- 
lutions (high shear modulus). We have not used more 
molecularly-motivated models; the GLAMM model is 
too complex to solve under shear banding conditions 
(IGraham et al. 1. 120031 iMilner. McLeish. and Likhtmanl . 
12001 ) , whil e the Rolie-Poly model has ze ro second nor- 
mal stress ( Likhtman and Graham! 120031) . 

Our model may thus explain some recently-reported 
phenomena: for example, redu cing the rim gap de l ays th e 
onset of meniscus instability ( Sui and McKennal 120071) ; 
a larger modulus G apparently increases the possibil- 
ity of instability (when comparing wormlike micelles to 
entangled polymer solutions); and there are variations 
in how far the applied shear rate progresse s alon g the 
stress plateau (jBerret. Porte, and Decruppel . 119971) . Fi- 
nally, our results allow one to rationalize the experiments 
on entangled polymer solutions that show a contentious 
mixture of shear banding and edge fracture: because of 
the relatively large second normal stresses of typical en- 
tangled polymer solutions, we naturally expect edge frac- 



ture to be associated with shear banding, which obviously 
complicates such measurements (as is apparent in the lit- 
erature). 
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Appendix A: Free Surface 

Here we outline the conditions necessary for a complete 
calculation of the free surface problem, which would de- 
termine the base state, including secondary flows, from 
which a true instability could be calculated. We consider 
steady creeping flow for which V • T = 0. Further, for 
planar Couette flow uniform in the x direction, stress 
gradients in x are zero. The free surface is subject to the 
boundary condition 



T ■ n — (-p a n) = ( 7s V ■ n)ii - V s7s , (Al) 



where 



and 



n 



V ■ n 



z-h'y 
VI + h' 2 

-h" 



(A2) 
(A3) 



is the curvature and h = h(y). We take the surface ten- 
sion to be uniform so V S 7 S = 0. 
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Eqs. (ED, (EU) and (pO| give 



x 

y 



-/^TQ^+T^ (A4) 
-h'T y y +^ y -JAtfM ■ n) , (A5) 
-h'T zy + %Y+ip* (A6) 



These three conditions determine the relation between 
the local curvature of the meniscus h'(y) and the stress 
components of the fluid evaluated at the meniscus. As 
h' — > the free surface becomes flat and we recover sim- 
ple shear conditions. However, generally h! can range 
from — oo to +00 as the shape of the meniscus changes, 
and this describes the limitation of our calculation. 

The remaining conditions are no flux through the free 
surface, v • n = 0; a stationary plate v(x, 0, z) = (0, 0, 0); 
and a moving plate v(x,W,z) — (V, 0, 0) . Far from the 
meniscus we require v(z — > —00) = (v x (y), 0, 0J,and the 
flow must reduce to simple shear flow given by uniform 
T yy and T xy , withT X2 = T zx = T yz = T zy = 0. From 
incomprcssibility V • v = 0, and the total stress is given 
by T = -pI + 2-qB + T,, where D = \ [Vv + (Vv) T ] and 
£ depends on the particular model. 

It is a challenging problem to find the flow and normal 
stresses at the free surface while also finding the shape 
of the fr ee surface, which is g e nerall y not circular. For 
example. iTanner and Keentokl ( 19831 ) assumed an initial 
semi-cylindrical crack in their calculation of edge frac- 

I 



ture, while iLoded ()l964fl showed that a spherical surface 
leads to a solution consistent with mechanical balance. 
Neither calculation will apply for a shear banded state 
where there must, by necessity, be different curvatures 
in the two bands due to different normal stresses. Addi- 
tionally, there is the complication of two (or more) bands 
with different shear rates in the bulk. 



Appendix B: Constitutive equations 

In this Appendix we present the homogeneous steady 
states of the two constitutive models we have used to 
fit the data in Section HVl We fit the total shear stress 



+ 777 to the measured flow curves of shear 
stress as a function of shear rate. Note that the flow 
curves include the stress plateau, while the constitutive 
curves are non-monotonic. We fix the position of the 
stress plateau based on the ex periments, rather than 
fittin g to a non-local model (|Lu. Olmsted, and Ball 
I2000D . This allows us to extract the second normal stress 
difference in the two shear bands, as the values ^2(71) 
and ^2(72) are determined by the shear rates 71 and 72 
in the coexisting shear bands, given by the intersection 
of the stress plateau with the two stable branches of the 
constitutive curve. 



1. Johnson-Segalman model 



In t he diffusive Johnson-Segalman (JS) model the viscoelastic stress X is assumed to obey ( Johnson and Segahnanl . 
119771 : lOlmsted. Radulescu. and Lul . l2000h 



(d t + v • V)S + (OS - Sfi) - a(DS + ED) + -S = 2-D + X>V 2 E, 

r t 



(Bl) 



where r is the linear relaxation time and a, which satisfies 
— 1 < a < 1, parametrizes slip of the polymer relative 
to the local flow field. The JS fluid has a Newtonian 
viscosity 77 due to the solvent and a polymer viscosity /x, 
which is related to the characteristic elastic modulus G 
by fj, = Gt. We define the viscosity ratio 



n 

e = — . 
A' 

Banding occurs for < a < 1 and < e < |. 
Variables are made non-dimensional according 



(B2) 



7 = T7, 



(B3) 



in terms of which the total stress, Eq. $Z§, is expressed 
as T = —pi + 2eD + X. We will use the same scaling for 



the Giesekus model. For planar Couette flow th e steady 
state s olution to the homogeneous JS equation is ([Larsonl . 



Il988t iBird. Armstrong, and Hassageii fl987t) 



Jxy 



J yy 



1 + (1 - a 2 ) 7 2 

(-i + «h 2 

l + (l-a 2 ) 7 2 ' 



(B4a) 
(B4b) 



in terms of which the second normal stress difference is 
given by N 2 = 



2. Giesekus Model 
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In the diffusive Giesekus model the viscoelastic stress £ is assumed to obey ( Giesekusl . ll982t iHelgeson et al. I. I2009h 

(B5) 



(d t + v- V)S + (OS - Ell) - (DE + SD) + -S = 2-D - -E 2 + £>V 2 £, 

r t jJL 



Here the non-monotonic behaviour comes not from slip, bu t from a non-linear relaxation term p a rametrized by a . 
The analytical solutions for planar Couette flow are given by (jBird. Armstrong, and Hassagerul987t lGiesekuall982h . 



where 



N 2 = 



A = 



1 - A 



l + (l-2a)A ' 



(1 - N 2 fj 
1 + (1 - 2a)A> 2 



1 + 16a(l - a) 7 2 - 1 



8a(l - a) 7 2 

A non-monotonic constitutive relation, and hence shear banding, is possible for a > \ and e < (2a — l) 2 /(2a) 



(B6a) 



(B6b) 



Appendix C: Three Band Configuration 
1. Meniscus profiles 

In this Appendix we compute the configurations for three bands, with a high shear band sandwiched between two 
low shear rate bands. Motivated by iBritton and Callaghanl 's experiments, we assume a central high shear rate band 
of width W2 and peripheral low shear rate bands w\ and ws, where w\ + w% + = W. An equivalent construction 
applies to a central low shear rate band. 

The height of the surface h(y) is given in the three regions by (Fig. 



h 1 (y)=R 1 



y_ 
Ri 



COS 0] 



(CI) 



' S?? 1 + i7 +C0S</)1 ) +( R i- R *)V 1 - + cos --Rising (C2) 



h 3 (y) = Ri 



,'y-W 

1 - I = COS 02 I - sin 2 



Ri 



H 



H=(R 1 - R 2 ) 



U>1 1^2 



COS 01 



and continuity of h' at y ~ w\ and y = wi + iu 2 requires 



IWl tU 2 W3 
— + — + — = - COS 01 - COS 02 • 
itl -«2 K3 



(C3) 



Continuity of h at y = w\ and y = W\ + W2 requires the difference in height H across the gap to be given by 

+ Ri (sin 02 - sin 0i ) , (C4) 



(C5) 



The normal stress balance conditions (Eq. [S]) require Rs = R\ . The width of the high shear rate band is given by 
W2 = (7app — 7i)/ (72 — 71 ) , while wi and W3 can, in principle, take any values as long as they satisfy W\ + W3 = 1 — Wi. 
Eqs. © and (JUS) lead to 



Ri 



1 



COS 0i — COS 4>2 — w 2 A 



and i?2 



■ COS 0i — COS 02 + (Wi + W^A 



(C6) 
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FIG. 14. (a) Profile of the meniscus for the three-band configuration for W2 = 0.45 and ib\ = 0.4, with contact angles cj>i — 140° 
and (f>2 = 130°, in which the outer bands have one shear rate 71 and the inner band has a different shear rate. Normal 
stress balance ensures that the two outer bands have equal radii of curvature. The width W2 of the inner band is given by 
u>2 = (iapp — 71)/ (72 —71)- (b) Energy as a function of ibi showing the separate contributions from surface tension (surface) 
and pressure (volume). The stable configuration occurs for wi ~ 0.52. 



2. Meniscus Integrity 

The conditions for integrity of the surface, h!(w\) < 00 and h'{w\ + w 2 ) < 00, lead to (from Eqs. IC6 

|(1 — Wi) COS 01 — Wi COS 02 — WiW2A\ < 1, 



(C7a) 



I [1 - (wi + W 2 )] COS (f>i - (wi + W2) COS 4>2 + U>2 [1 - (wi + w 2 )] A\ < 1. 



(C7b) 



These conditions apply at both interfaces of the central 
high shear rate band, and must be satisfied simultane- 
ously for the meniscus to maintain integrity. The contact 
angles <f>i, (j) 2 and the distortion parameter A are prop- 
erties of the fluid, while the width of the high shear rate 
band w 2 is determined by the applied shear rate. 

Although a symmetric solution w\ — W3 is appealing, 
we will consider the more general case where the two 
low shear rate bands need not be the same size. We 
will see that a simple stability analysis can lead to sym- 
metry breaking to a non-symmetric band configuration. 
As has already been mentioned, stability should strictly 
be determined by dynamical considerations. However, 
since we are only treating the mechanics of the meniscus, 
we will study the three band configuration by construct- 



I 

ing an energy function, which roughly corresponds to the 
work needed to establish an interface, and whose gradi- 
ent specifies a generalized force that should vanish for 
mechanical stability. 

The mechanical energy of the surface, per unit length 
in the flow direction, is given by 



£(wi,W2, A,4>-l,4>2) = 7 S / vi + h' 2 dy - Ap / h dy 
Jo Jo 

(C8) 

which comprises the surface tension and the work done 
in deforming against the pressure differences across the 
interface. We take £ = at the onset of banding. For a 
given width of the high shear rate band W2 we determine 
the position of the band by minimising £ with respect to 
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w\. We will only consider the equal contact angle case, 
4>\ = 4>2 = 4>- Stable configurations are defined as those 
for which the energy function £ is a local minimum with 
respect to changing the sizes w\ and W2 of the two outer 
shear bands at fixed W2- 



3. Different three-band configurations 



Three configurations of the three band model minimize 
the energy, as demonstrated in Figs. 1151 



(i) central high shear rate band (w\ =1(73); 

(ii) off center high shear rate band (w\ 7^ W3); 

(iii) collapsed low shear rate band (wi = 0, or W3 = 0), 
corresponding to two shear bands. 

We can show that, for equal angles, there will be an ex- 
ternum at W\ = 2(1 — ^2)- The high shear rate band 

remains central so long as — — % 



> 0. Equality de- 
|(i-w 2 )_ _ 

fines W2, so that the condition for stability is W2 < W2 
where u>2 satisfies 



1 - 



27.', 



A w*2 Ri 



2 cos <b W R 2 



Ri 
R> 



(0 



0. 



(C9) 



1 — ( cos <ft 
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together with the integrity conditions of equation (|C7|) . 
Further, for given distortion parameter A and high shear 
rate band width W2 , a three band configuration may sat- 
isfy the integrity conditions of Eq. (|C7|) . At a slightly 
greater value of W2 (or applied shear rate) the three band 
configuration becomes unstable and the fluid collapses 
into a two-band states, with integrity governed by equa- 
tion (fT3)) . 

Fig. Q2] shows the energy curves and some meniscus 
shapes for the distortion parameter A = —3.5 and equal 
contact angles <fr = 130°, as a function of increasing 
the width of the central high shear rate band W2- For 
^2/(1 — W2) = 0.5 the low shear material splits into two 
equal size bands surrounding the high shear rate band. 
Upon increasing 1I12 the symmetric three band configu- 
ration becomes unstable to an asymmetric configuration 
(for u>2 > 0.343), and finally unstable to a two band con- 
figuration when the energy minima lie at either w± = 
or iii = 1 — u>2- 

Fig.lTSlshows all loci of states whose meniscus integrity 
we can estimate using this method for A = —3.5. This 
figure shows the regions, in Fig. [SJ in which two bands 
fracture at different ends of the stress plateau depend- 
ing on whether the contact angle is closer to wetting or 
non-wetting angles. In addition, there are accompanying 
regions in which the inner of three bands is either sym- 
metric or non-symmetric. In this particular example, the 
larger contact angles (<fi > 90°) correspond to a central 
high shear rate band, while the smaller contact angles 
(4> < 90°) correspond to a central low shear rate band. 
There are four distinct classes of regions of for which 
the meniscus retains integrity: (1) two bands, (2) three 
bands with a central symmetric band, (3) either three 
symmetric bands or two bands, or (4) either three asym- 



metrically distributed bands or two bands. In the last 
two cases where two band configurations allow meniscus 
integrity we expect factors such as flow history, flow ge- 
ometry (e.g. cylindrical Couette vs cone and plate) and 
boundary conditions to determine the selected configu- 
r ation. 

iBritton and Callaghanl (|l999T l found the central high 
shear rate band to be stable and to increase in width as 
the applied shear rate increased (Fig. [17]). At 7 = 7 s -1 
the high shear rate band is stable in a central position 
consistent with Fig. [T5b . while at 7 ~ 10.7 s" 1 the high 
shear rate band is no longer central but in an off-center 
stable position consistent with Fig. [15b . This would be 
consistent with the state diagram in Fig. [1^1 for <f> > 120°. 



4. Summary 

In summary, we find the following results for three 
bands: 

1. Three bands are stable in a subspace of the the 
parameter space, for a range of contact angles. This 
stability region starts with an infinitesimally small 
central band, and continues until the central band 
becomes too large, after which only a conventional 
two band state is possible. 

2. The stable central band is the high (low) shear rate 
branch when the contact angle is closer to 180° (0°). 

3. The central of the three bands is symmetric for 
smaller central bands, but is destabilized towards 
an off-center configuration for wider central bands, 
for a range of contact angles. 
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FIG. 15. Total energy curves (a) and mensci shapes (b,c,d) for distortion parameter A = —3.5 and contact angle (f> = 130°. 
For widths W2 < 0.343 of the high shear rate bands (a), the central high shear rate band is stable. For u)2 > 0.343 the high 
shear rate band is no longer stable at the central position but stabilises such that there is a wide low shear rate band and a 
narrow low shear rate band close to the plate (c,d). For ui2 > 0.579 the fluid collapses into the two band configuration. The 
meniscus profiles are shown for W2 — 0.2, 0.4, 0.5, and the complete set of states for all contact angles is shown in Fig. 1161 



4. In some regions of the phase space two-band 
and three-band states are simultaneously sta- 
ble; in other regions the two band state is un- 
stable, in principle towards three band states. 
This may explain some of the experiments of 
iBritton and Callaghanl (|l997t) . in which three band 
states are found. 

5. We are not able to assess the relative stability of 
two or three band states, and we speculate that 
this is determined by factors such as initial condi- 
tions, boundary conditions, and the degree of cur- 
vature in the flow. For example, the highly curved 
flow of a wide gap Couette rheometer may have 
a strong preference for the two band state, as it 
also leads to a preference for inducing the high 
shear rate phase next to the inner cylinder where 
the shear stress is highest dRadulescu and Olmstedl . 
l200fj lAdams. Fielding, and Olmstedl . 120081 ). 
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FIG. 16. Integrity diagram for A = —3.5, show regions in which two or three bands allow meniscus integrity as a function of 
contact angle and the fraction of material W2 in the middle band. For <j> > 90° (<f) < 90°) the central band adopts the high 
(low) shear rate. 




FI G. 17. NMR images of a worm -like micellar solution undergoing shear banding in a cone-and-plate rheometer (image from Fig. 
lO lBritton and Callaghanl l|l999h ). At an applied shear rate of j = 7 s _1 the high shear rate band, appearing as white/black, 
is stable and central. This would correspond to Figure (|15b~l. At the higher applied shear rate 7 = 10.7 s~ , the high shear 
rate band is stable but no longer central but positioned close to the moving cone. This corresponds to Figure (|15b ). At even 
higher shear rates the material was expelled from the cell. 
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